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We present the isospin- and angular-momentum-projected nuclear density functional 
theory (DFT) and its applications to the isospin-brcaking corrections to the superallowed 
/3-dccay rates in the vicinity of the N = Z line. A preliminary value obtained for the 
Cabbibo-Kobayashi-Maskawa matrix element, \V u d\ = 0.97463(24), agrees well with the 
recent estimate by Towner and Hardy [Phys. Rev. C77, 025501 (2008)] . We also discuss 
new opportunities to study the symmetry energy by using the isospin-projected DFT. 



1. Introduction 

The isospin-symmetry violation in atomic nuclei is predominantly due to the 
Coulomb interaction that exerts long-range polarizations on neutron and proton 
states. To consistently take into account this polarization, one needs to employ huge 
configuration spaces. For that reason, an accurate description of isospin impurities 
in atomic nuclei, which is strongly motivated by the recent high-precision measure- 
ments of the + — > + Fermi superallowed /3-decay rates, is difficult to be obta ined 
in shell-model approaches, and specific approximate methods are required!!^ 

The long-range polarization effects can be included within the self-consistent 
mean-field (MF) or DFT approaches, which are practically the only microscopic 
frameworks available for heavy, open-shell nuclei with many valence particles. These 
approaches, however, apart from the physical contribution to the isospin mixing, 
mostly caused by the Coulomb field and, to a much lesser extent, by isospin- 
non-invariant components of the nucleon-nucleon force, also intro duce th e spurious 
isospin mixing due to the spontaneous isospin-symmetry breakingPEEl 

Hereby, we present results on the isospin mixing and isospin symmetry- 
breaking corrections to the superallowed Fermi /3-decay obtained by using the 
newly deve loped i sospin- and angular-momentum-projected DFT approach with- 
out pairingPEEMl xhe model employs symmetry-restoration techniques to remove 
the spurious isospin components and restore angular momentum symmetry, and 
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takes advantage of the natural ability of MF to describe self-consistently the sub- 
tle balance between the Coulomb force making proton and neutron wave functions 
different and the isoscalar part of the strong interaction producing the opposite 
effect. 

The paper is organized as follows. In Sec. [2] we describe the main theoretical 
building blocks of the isospin- and angular-momentum-projected DFT. Section [3] 
presents some preliminary applications of the formalism to the isospin symmetry- 
breaking corrections to the Fermi superallowed /3-decay matrix elements, whereas 
Sec. H] discusses applications of the isospin-projected DFT to nuclear symmetry 
energy. The summary is contained in Sec. [5] 

2. The projected DFT framework 

The building block of the isospin-projected DFT is the Slater determinant, $), 
represen ting the self-consistent Skyrme-HF solution provided by the HF solver 
HFODDP^ Self-consistency ensures that the balance between the long-range 
Coulomb force and short-range strong interaction, represented in our model by 
the Skyrme energy density functional (EDF), are properly taken into account. The 
unphysical isospin mixing is taken care of by the rediagonalizati on o f the entire 
Hamiltonian in the good isospin basis, \T,T Z ), as described in RefsP^I This yields 
the eigenstates: 

\n,T z )= <4,tJT,T z ) (1) 

T>|T 2 | 

numbered by an index n. The so-called isospin-mixing coefficients (or, equivalently, 
isospin impurities) are defined for the n— th eigenstate as 

a C = 1 ~ l a T,Tjmax> ( 2 ) 

where |oij. T ^|^ lax stands for the dominant amplitude in the wave function \n,T z ). 

Within the isospin- and angular-momentum-projected DFT, we use the nor- 
malized basis of states \I,M,K;T,T Z ) having both good angular momentum and 
good isospinfm Here, M and K denote the angular-momentum components along 
the laboratory and intrinsic z-axes, respectively. The K quantum number is not 
conserved. In order to avoid problems with overcompleteness of the basis, the K- 
mixing is performed by rediagonalizing the Hamiltonian in the so-called collective 
space, spanne d for each I and T by the natural states, \IM;TT z )^ l \ as described 
m Refs.™ Such a rediagonalization yields the eigenstates: 

\n;IM;T z )= £ a^ T \IM;TT z f\ (3) 

i,T>\T,\ 

which are labeled by the index n and by the conserved quantum numbers /, M, 
and T z = (N - Z)/2 [compare Eq. JTJ]. 

The isospin projection does not produce singularities in energy kernels; hence, 
it can be safely used with all commonly used EDFs. Coupling the isospin and 
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Fig. 1. The absolute values of the norm kernels, |A/"(/?t; a , 7) = | {^l -R(/3t)-R(c*, /3, 7)1*) I, for 
a state in 14 N calculated with the SLy4 EDF, plotted versus the rotation angle in the isospace 
The solid curve, exhibiting the single singularity at /3y = 7r, corresponds to the pure isospin- 
projected DFT theory, which is regular for all Skyrme-type functionals.El The dotted lines corre- 
spond to two fixed sets of the Euler angles in space, with a = 7 ~ 0.314, and /3 « 0.229 (left curve) 
and j3 ~ 1.414 (right curve). The poles that appear inside the integration region, < /?t < 
give rise to singularities in the projected DFT approach. 



angular-momentum projections, however, leads to singularities in both the norm 
(see Fig. Q]) and energy kernels. This fact narrows the applicability of the model 
to Hamiltonian-driven EDFs which, for Skyrme-type functionals, leaves only one 

[To] 

option: the SV parametrizationM^ The alternative would be to use an appropriate 
regularization scheme, which is currently under development ! ^ ^ ^ * 



3. Isospin-mixing and isospin-breaking corrections to superallowed 
(3- decay 

Evaluation of ac is a prerequisite to calculate isospin corrections to reaction and 
decay rates. As is well knownpl isospin impurities are the largest in N — Z nuclei, 
increase along the N — Z line with increasing proton number, and are strongly 
quenched with increasing \T Z \ — \N — Z\/2. Such characteristics were also early 
estimated based on the perturbation theorjOSI or hydrodynamical modelP^ Quan- 
titatively, after getting rid of the spurious mixing, which lowers the true ac by as 
much as 30% P the isospin impurity increases from a fraction of a percent in very 
light N = Z nuclei to -0.9% in 40 Ca, and -6.0% in 100 Sn, as shown in Fig. El In the 
particular case of 80 Zr, the calculated impurity of 4.4% agrees well with the empir- 
ical value deduced from the giant dipole resonance 7-decay studies.^ This makes 
us believe that our model is indeed capable of capturing essential physics associated 
with the isospin mixing. Unfortunately, current experimental errors are too large 
to discriminate between different parametrizations of the Skyrme functional. The 



26, 2010 1:32 WSPC/INSTRUCTION FILE kazilO-08 



4 W. Satula, J. Dobaczewski, W. Nazarewicz, M. Borucki, M. Rafalski 



i_ ., ( y = 2.6787 - 0.055649X R= 0.77657 ) 



i£ 6 
4 




( y = 24.1 93 - 0.54926X R= 0.91 273] 



31.5 32.0 32.5 33.0 33.5 34.0 34.5 
doorway state energy [MeV] 



Fig. 2. Isospin impurities in the ground states of 40 Ca (upper panel) and 100 Sn (lower panel), 
plotted as fun ctio ns of the excitation energy of the doorway state for a set of commonly used 
Skyrme EDFs.DSl Results of the linear fits and the corresponding regression coefficients, R, are 
also shown. 



variations between EDFs in Fig. [2] result in ~10% uncertainty in calculated values 
of ac- 

The magnitude of theoretical ac is quite well correlated with the excitation 
energy, £t=i> of the T = 1 doorway state, see Fig. [2] However, in order to make 
a precise determination of Et=i, spectroscopic quality EDFs are needed, and this 
is not yet the casePH This explains why the values of ac do not correlate well 
with basic EDF characteristics, including the isovector and isoscalar effective mass, 
symmetry energy, binding energy per particle, and incompressibility (see discussion 
in RefP). 

Increasing demand on precise values of isospin impurities has been sti mulat ed 
by the recent high-precision measurements of superallowed /3-decay rates.^IH A 
reliable determination of the corresponding isospin-breaking correction, 8c, requires 
the isospin- and angular-momentum-projected DFtM This correction is obtained 
by calculating the + — > + Fermi matrix element of the isospin raising/lowering 
operator T± between the ground state (g.s.) of the even-even nucleus \I = Q, T w 
1,T Z = ±1) and its isospin-analogue partner in the N = Z odd-odd nucleus, \I — 
0,T«1,T Z = 0): 

|(/ = 0,Tw l,T z =±l|f±|J = 0,Twl,T« =0)| 2 =2(1 -8c). (4) 

To determine the \I — 0,T m 1,T Z = 0} state in the odd-odd N = Z nucleus, 
we first compute the so-called antialigned g.s. configuration, |j>®7r) (or |^®7f)), by 
placing the odd neutron and the odd proton in the lowest available time-reversed 
(or signature-reversed) HF orbits. Then, to correct for the fact that the antialigned 
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Fig. 3. Values of \V ul i\ deduced from the superallowed /3- decav (full circles) for three different 
sets of the 5q corrections calculated in: RefPl (a); Ref|22l with NL3 and DD-ME2 Lagrangians 
(b); and in the pr esen t work (c). Triangles mark values of \V u d\ obtained from the pion-decajED 
and neutron-decaj!24J studies, respect ivel y. The open circle shows the value deduced from the 
/3-transitions in T = 1/2 mirror nucleiPS 



configurations manifestly break the isospin symmetrypl that is, \i? n) w (|T = 
0) + \T = 1}), we apply the isospin and angular-momentum projections to create 
the basis \I, M, K, T, T z = 0), in which the total Hamiltonian is rediagonalized (see 
Sec. [2]). A similar scheme is used to compute the \I = 0, T w 1, T z = ±1} states in 
the even-even nuclei. 

Our studies indicate^ that to obtain a fair estimate of 5c for A < 40 and A > 40 
nuclei, one needs to use large harmonic oscillator bases consisting of at least TV = 10 
and 12 full shells, respectively Even then, the results are subject to systematic er- 
rors due to the basis cut-off, which can be estimated to be ^10%. Despite the fact 
that not all TV — 12 calculations in heavy (A > 40) nuclei have yet been completed, 
and that owing to the shape-coexistence effects, there are still some ambiguities 
concerning the global minima, our preliminary results point to encouraging conclu- 
sions. Namely, the mean value of the structure-independent statistical-rate function 
.Fip^ obtained for 12 out of 13 transitions known empirically with high precision 
(excluding the 38 K— >- 38 Ar case), equals J-t — 3069.4(10), which gives the value of 
the CKM matrix element equal to \V u d\ — 0.97463(24). These values match well 
those obtained by Towner and Hardy in their recent compilation^ (see Fig. [3]) . Be- 
cause of a poor spectroscopic quality of the SV parameterization, the confidence 
levePZl of our results is poor. Nevertheless, it should be stressed that our method 
is quantum-mechanically consistent (see discussion in Refsl ^^ j) and contains no 
adjustable free parameters. 

4. Symmetry energy 

The spontaneous violation of isospin symmetry in all but isoscalar MF configura- 
tions ol N = Z nuclei offers a way to study the nuclear symmetry energy. The idea, 
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Fig. 4. Top: schematic illustration of the isospin-symmetry-breaking mechanism in MF of odd-odd 
N = Z nuclei. Bottom: i?iym' m odd-odd N = Z nuclei calculated with SLy4, SV, SLy4^, and 



which is schematically sketched in the upper portion of Fig. @J invokes the mixed- 
symmetry antialigned \v ® 7r) (or \v (g> 7f}) configuration in an odd-odd N = Z 
nucleus. By applying the isospin projection to the HF state \v <8> 7r), one decom- 
poses it into the isoscalar T = and isovector T = 1 parts. As argued below, the 
magnitude of the splitting, I^lym , depends on the isovector channel of a given EDF, 
i.e., its symmetry energy. 

For the Skyrme- ty pe EPFs, the symmetry energy in the nuclear matter limit 
can be decomposed as!^ 



x sym 



1 

t: £ fg 



3tt 2 \ 2/3 



= a (kin) 

sym 



,(int) 



(5) 



The first term in Eq. ([S]) is associated with the isoscalar part of the nucleon-nucleon 
interaction and primarily depends on the mean single-particle level spacing at the 
Fermi energy. This term is scaled by the inverse isoscalar effective mass. The second 
(interaction) term, is related to the isovector part of the Skyrme-EDF: 8% t= i = 
C{p\ + C\p\T\ (for definitions, see RefP 

and references quoted therein). 
The value of i?iym^ appears to be mainly sensitive to the interaction term, which 
is illustrated in Fig. 01 Indeed, despite the fact that SLy4 and SV EDFs have similar 
values of a sym (equal to 32MeV and 32.8 MeV, respectively), the corresponding 



energy splittings -E S ym^ differ substantially. The reduced values of |£7syrn| m SV are 
due to its small value of esi^ = 1.4MeV|fJ which is an order of magnitude smaller 



(int). 



a This small value shows how unphysical are the consequences of the saturation mechanism built 
into SV through the strong momentum dependence and results in an unphysically low isoscalar 
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than the corresponding SLy4 value: aiym = 14.4 MeV. 

An interesting aspect of our analysis of i?iym' relates to its dependence on 
the time-odd terms, which are poorly constrained for Skyrme EDFs. To quan- 
tify this dependence, we have performed calculations by using the SLy4^ and 
SkM^ functiona ls, wh ich have the spin coupling constants adjusted to the Lan- 
dau parametersP^Hl] These EDFs have different values of a sym but the same 
Osym = 14.4 MeV. The similarity of the calculated energy splittings shown in Fig.rj] 
confirms that this quantity primarily depends on the isovector terms of the func- 
tional. Moreover, its significant dependence on the time-odd terms opens up new 
options for adjusting the corresponding coupling constants to experimental data. 
This will certainly require the simultaneous restoration of isospin and angular- 
momentum symmetries, as presented in this study. 

5. Summary 

In summary, the isospin- and angular-momentum-projected DFT calculations have 
been performed to estimate the isospin-breaking corrections to + — > + Fermi 
superallowed /3-decays. Preliminary results for the average value of the nucleus- 
independent Ft = 3069.4(10) and the amplitude \V ud \ = 0.97463(24) were found to 
be consistent with the recent estimates by Towner and Hardypl notwithstanding a 
low spectroscopic quality of the Skyrme EDF SV used. 

Applicability of the isospin-projected DFT to analyze the nuclear symmetry 
energy has also been discussed. It has been demonstrated that the isospin projection 
offers a rather unique opportunity to study the interaction part of the symmetry 
energy in the odd-odd N = Z nuclei and that this quantity is influenced by time-odd 
fields of the energy density functional. 
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